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Der Aufbau und erste Versuchsergebnisse von neuartigen, ganz aus 
Kunststoff gefertigten, flexiblen Solar-Luftkollektoren mit porösem Absorber 
sowie das hierzu entwickelte Meßverfahren werden beschrieben. Typische 
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Angewandte Systemanalyse 












































Anwendungsmöglichkeiten für diese Kunststoffkollektoren liegen im Agrarbereich, Aan begt design of solar air heaters of porous type made completely 3.1 € for framing the collector 
wie die Trocknung von Ernteprodukten, jedoch auch Brauchwassererwarmung und p astic materials has been forwarded for typical applications like I roms 4 
Nr. 26 Raumheizung gehören dazu. crop drying, domestic water and space heating etc. An experimental test 3,3 Air flow measurement 9 
j facility brought up to study the thermal performance of these collectors Wr MD LL d 
i i Für unsere Versuche wurden drei verschiedene Kollektoren dieses Typs has been described. 3.5 E sna BR 
Plastic Solar Air Heaters angefertigt, von denen zwei eine Lange von etwa 10 m und eine Breite von Three collect i ٧ Wi à ape a sasurements 15 
1 m haben, sich jedoch im Kunststoffgewebe des Absorbers, das von vers.hiedenen cated for Pond ion ors, each having a width of about one meter, were fabri- D e وېي‎ qi سر‎ ۳ 

of a Novel Design Firmen stammt, unterscheiden. Der dritte Kollektor ist 20 m lana und vurde esting purposes. Two of these collectors have identical length "m raai " 


of approximately ten meters but the absorber, which is a | 

obtained from two differrent firms. The third collector Mr ومو‎ 
long and was essentially made to study the effect of length. A detailed 
study was undertaken to determine the variation in the thermal performance 
of these collectors with typical parameters like quantity of the flowing. 


air, solar insolation, absorber type and back insulation, etc. For 
identical ambient conditions and for a fixed air flow, the performance 


of the collectors is found to be dependent very sersitively on the ki 
textile absorber and the insulation at the back. ru 


vermessen, um den Effekt unterschiedlicher Kollektorlängen zu untersuchen. 














- Testing and Performance - 








Der thermische Wirkungsgrad der Kol'^'toren wurde in Abhängigkeit von ٩ —— — de 
den typischen Parametern, Luftmenge, solare Einstrahlung, dem Absorbegewebe- 
typs sowie der Dicke der rückseitigen Isolierung gemessen. Um ein 5۹ 
Verständnis der Luftströmung und Wärmeübertragung im Kollektor bei den 
verschiedenen Absorbegeweben zu gewinnen, wurden die Luftdruck- und Tempera- 
turverteilung über die Länge oberhalb und unterhalb des Absorbertuches 
gemessen. 
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Für identische Umgebungsbedingungen und identische Luftmenge zeigt sich F i 6.1 10 m Collector a 49 

der Wirkungsgrad stark abhängig von der Art des Absorbegewebes und der rück- , For the better 10 m long collector a high thermal performance is 6.2 10 m Collector b 27 

seitigen Isolierung des Kollektors. Die besten Ergebnisse lassen sich mit achieved if one keeps a 6 cm thick polyethylene insulation at the back of 6.3 20 m Collector 30 

einem der 10 m Kollektoren bei einer 6 cm starken rückseitigen Isolierung aus the collector. For an incident solar radiation of 687 w/m" at the collector 6.4 Comparison of the Collectors 33 

wasserabweisende Polyäthylenenschaum erzielen. Bei einer solaren Einstrah- surface, one is able to achieve a temperature rise of 16.6° C in the air 6.5 Effect of Back Insulation 35 

lung von 687 W/m? in der Kollektorebene ergibt sich bei einem Volumenstrom flowing through the solar collector at a rate of 800 m'/h, thus yielding 6.6 Presente Step la Che یي وغ‎ di 

| - 8 i an efficiency of nearly 71 % and useful 1 6.7 Effect of Relative Humidity on the Thermal Efficiency 38 

von 800 m'/h Luft durch den Kollektor eine Temperaturerhöhung von 16,6 C ermal power of about 4.5 kW. 240200 وښ‎ Seétel-dBiles Semer Ae n 
und damit ein Wirkungsgrad von etwa 71 % oder eine nutzbare thermische 6.8 Validity of Ho q 















Further it is seen that the linear approximation for the Hottel- 
Bliss equation leads to erroneous estimations for the coliector parameters 
like the total heat loss coefficient U » plate efficiency factor F' and the 
effective (qt). product, when the absbrber is porous. For a dense absorber 
however, such ah approximation yields, as usual, correct numerica] values i 
for the characteristic parameters of the collectors. 


Leistung von 4,5 KW. 























Die theoretische Betrachtung des Kollektors zeigt, daß die in der 1 — 
Temperatur lineare Hottel-Whillier-Bliss Gleichung zu falschen Werten für 
die Kollektorgrößen Wärmeverlustkoeffizient U, , effektiver Transmissions- 
Absorptionsfaktor (et) und Kol lektorwirkungSgradfaktor F' führt, wenn es 
sich um stark poröse AbSorberstrukturen handelt. Für das dichtere Absorber- 
gewebe führt die lineare Näherung, wie gewöhnlich, zu korrekten Werten für 
die charakteristische Größen des Kollektors. 
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1 INTRODUCTION 


The potential applications of solar air-heaters are domestic 
hot water and space heating, drying of agricultural products and lumber 
under controlled conditions, heating of greenhouses etc. In contrast to 
vast studies on water heating collectors, the studies on the design and 
performance of air heaters have been rather limited. The difficulties of- 
ten encountered with air heating collectors are the limitations of heat trans- 
fer from the absorber plate to the working fluid and the need to handle 
large volumes of air than liquids due to the low density of air as the work- 
ing substance. Several designs of solar air heaters have been proposed 
for improving their thermal performance and some of them have yielded good 
results /1, 2/. 


ne designs of solar air heaters forwarded so far can basically 
be classified into two categcries. The first type has a nonporous absorber 
in which the air stream does not flow through the absorber plate. Air may 
flow above and/or behind the absorber plate. The second type has a porous 
absorber that includes slit and expanded metal, transpired honey-comb and 
overlapped glass plate absorber. The utility of the later type over the former 
one consists in providing a lower pressure drop across the collector and 
also reducing thermal losses to the environment because the solar radiation 
in this case penetrates to greater depths and is absorbed gradually depend- 
ing upon the matrix density /1/. The studies on nonporous air heaters have 
been made by Whillier /3/, Close /4/, Khanna and Singh /5/, and Gupta and 
Garg /6/, Bansal and Kaushik /7/ and Parker /8, 9/ etc. Porous matrix type 
air heaters have been extensively analysed and tested by Chiou et al /10/, 
Beckmann /11/, Hamit and Beckmann /12/, Selcuk /13, 14/, Akyurt et al /15/, 
Lansinget al /16/, Singh et al. /17/ and Sodha et al /18/. 
In most of these studies the material used for the absorber plate was that 
of metal and for glazing one used either glass or plastic sheets. The ~ost 
of the collector therefore did not come out to be much cheaper than those 
of water heating collectors. 


One of the ways to make the solar heating economically attrac- 
tive is to use the plastic materials in the manufacture of various components. 
A review listing various plastic materials suitable for various components 
of a solar energy collector has been given by Blaga /19/ and Boeckmann 
/20/. The procedure of making thin films of the plastic on the absorber and 
over the cover has been described by Andrews and Wilhelm /21/. The solar 
air heating collectors made completely from the plastic material include 
those of Wilhelm /21/, Weichman and Hughes /22/, Sarazin and Olson /23/, 

Chio /24/, Erb /25/, Hoynowski /26/, Kay /27/, Boettcher and Delyannis /28/ 
and Schulz /29/. The plastic collectors constructed from the plastic materials 
so far have been mainly of nonporous type which can abviously reach a limited 
efficiency and have the disadvantage of providing high pressure drop across 
the collector. 


In this communication we present a detailed study of porous 
type of air heating solar collectors made from piastic material and which 
are capable of handling large volume of air. The design of the collector 
has been described in Section 2. Three different collectors, two of 10 m 
length and another of 20 m length have been fabricated by Dynamit Nobel. 


m 





t Dynamit Nobel, 521 Troisdorf, West Germany 
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The detailed dimensions along with typical physical parameters 
of the three collectors are described in table | of section 2. To test the 
collectors an outdoor facility was set up at the terrace of the institute 
(STE) and the measurements of the reievant quantities were done on a MADAS 
/ 33. 37/ data processing system. The details of the test facility and ex- 
perimental set up are described in section 3. 


Section 4 deals with the relevant performance equations and in 
section 5 we present the results of measurements and calculations. 


DESIGN OF THE COLLECTORS 


The present collector belongs to the category of porous absor- 
bers and is fabricated from flexible plastic sheets. As an absorber there 
is a porous black textile of polyester (100 Z) and it is attached to two 
transparent sheet of polyviny! chloride (PVC) at the edges and is covered 
by them from both the sides (Fig. 1). The edges of the collector have ring 
hoies at a distance of 0.50 m from each other for fixing the collector on 
the stand (described in section 3). The transparent PVC sheet which acts as 
cover is 0.6 mm thick and has about 83 Z transmittance for global solar 
radiation and the material is stabilized against ultraviclet radiation. 

The measured transmittance* of the present PVC sheet in the spectral range 
9-20 mm is very low (~ 6 2) (Fig. 2). Thus for absorber temperatures in the 
range 40 - 100 9C this PVC sheet acts as an excellent cover. The lower 

PVC sheet is provided to double the life time of the collector; when the 
top cover becomes opaque to solar radiation due to natural effects then 

the collector is simply inverted so that the bottom PVC sheet now acts as 

a cover. 


The porous textile in the collector should be so selected as 
to have maximum absorption of solar radiation but allowing a good flow of 
air through it. Since the air flows over the entire absorber area the pres- 
sure drop across the collector in such constructions is much lower than in 
the case of nonporous absorbers. Simultaneously, since the radiation is ab- 
sorbed by different layers over the thickness of the absorber, the upper 
layer is not at the maximum temperature resulting in reduced thermal losses. 
Presently we got three collectors fabricated. Two of them are of approxi- 
mately 10 m? area each and the other is of about 20 m? area. The width of 
the black textile in each case is about 1 m, while the top and bottom PVC 
covers are made slightly bigger. With this construction the absorber lies 
flat after fixing while the air is contained between the top PVC sheet and 
textile respectively, before and after getting heated. Our two 10 m area 
collectors have two different textiles found most suitable from preliminary 
measurements. A twenty times enlarged surface photograph of the two textiles 
taken from an electron microscope is given in Figures 3a and 3b respective- 
ly. One 10 m2 area collector (called 10 m collector a later) has the textile 
of Fig. 3a as the porous absorber. This textile was found to transmit 
about 8 % of the solar radiation. The absorptance of the absorber is 


therefore aro mote ty 88 % assuming 4 % reflection losses from its surface. 


Another 10 m^ area collector (called 10 m collector b 





+ 
Transmittance measurements were performed at the Indian Institute of 


Technology, Delhi, India 


textile absorber 


PVC 


Fig.1: Diametric section of the air collector 


Fig.2: Measured transmittance of PVC-cover in the wavelength region of 
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2.5 - 25 pun (cover thickness O.6 mm) 
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Fig. 3a: Surface structure of the textile absorber of 10 m collector a 


(20 times magnifiod electron microscope photograph) 


Fig. 3b: Surface structure of the textile absorber of 10 m collector b and 20 m 
collector (20 times maonified electron microscope photograph) 


(e.m. photographs of the textile surface were made at IGV/KFA-Jülich) 
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later) and 20 m? area collector (called 2 
tile of Fig. 3b as the absorber. This text 
78 % of the solar radiation. The details 

their geometry and relevant physical Para 


Table 1: Dimensions and physical properti: 








Items Om collector : 
length of the collector 8.85m 
width of the collector 1.03m 
absorbing area AE 


cross-section of the collector 


| 0.061m? 
at the inlet and outlet of th 
air 
material of the cover PVC 
thickness of the cover 0.6mm 


transmittancet of the PVC- 
cover in the solar spectrum 0.83 
range 


transmittance in thc wave- 


length range 9-20ym 0.06 

material of the absorber 100% Polyester 

tissue textile No.: 
12503* 


thickness of absorber tissue |~2mm 


transmittance of the textile 


absorber for solar rediation | 0.08 











** Firm. Säckinger Feinwirkerei,7880 Säckinger 


O m collector later) have 
ile (Fig. 3b) is found to 
of each of the collectors 
meters are given in table 


25 of the plastic collectors 


| |10m collector b 


8.80m 17.60m 


1. 15m 1.13m 
10. 12m? 19.89m? 


0.061m2 0.061m2 


PVC PVC 
0.6mm 0.6mm 


0.83 0.83 


0.06 0.06 


100% Polyester |100% Polyester 


textile No.: Textile No.: 
2262/180 ** 2262/180** 


~ 2mm ~~ 2mm 


0.22 0.22 


/Rhein,FRG 





EXPERIMENTAL SET UP AND MEASURING INSTRUMENTS 


The test facility for testing the solar air collector is schema- 
tically shown in Fig. 4. It consists of the following set up. 


pyranometer 


aperture 


blower 


Fig.4: Schematic of Test Facility 


Stand for framing the collector 
i Si 


The collector being flexible has to be put very properly and 
stretched in such a position so as to receive maximum solar radiation. The 
frames on which the collector was to be fixed were 5 meters long and 1.50 m 
wide. The frames were fixed on 5 tables oriented towards south with the 
possibility of changing the angle which was kept 459 for the present study.** 
They were covered with wood and polyethylene insulation of 3 cm thickness 


initially. The collector was fixed over the insulation. A photograph of the 
completed test facility is given in Fig. 5. 


Jed Blower 


The blower used for letting the air flow inside the collector E = 


is a Lau forward curved blower /30/ from the firm Fischbach* of type D570/ 
E280. The details of the blower are as follows: 


IN‏ سس 
+ 


Fischbach, D-5908 Neunkirchen, West Germany 
+ 


Latitude of Jülich is ¢ = 500 55' 33" north 
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VIVAL DUMEPU U ANDAHDS 
STANDARD REFERENCE MATERIAL 1010a 
(ANSI and ISO TEST CHART No. 2) 


Response time | second (l/e signal) 


DIN* 19205 /31/. The introduction of a multinozzle plug between the manome- 


ter and the nozzle enables to measure the static pressure also. Cosine | per cent from normalization O - 


zenith angle 


The amount of air m flowing through the collector is calculated 3 per cent 70 - 80" zenith angle 


by using the following equation 


Ori : : : 
Absclute electrical 0.25 kW rientation no effect on instrument performanc 
power h=8.€. A . V2 AP.e | M : : : d die. دند‎ 

Electrical potential variation 30 - 220 V Y w d ۶ : | echanıcal vibration tested up to 20 g's without damage 
E for f ti dl ZEN : Calibration integrating hemisphere (approx. | 

x. Temp. for free operation | where the flow rate number d and the expansion factor £ both are functions | cm ^ min ', ambient temperature +2 

Max; 1 e blowi TU of the Reynold number and the dimensions of the measuring track. For the calibration reference Eppley prima 

—— — سم‎ 7 "- RW "e: em | aperture supplied to us and for our measuring track m is simply calculated dard group of Angstrom pyrheliomet 

| by the formula dag reproducing the International Pyrh 


max 
metric Scale 


Readout Potentiometric in preference to mi 


Air flow measurement ammeter devices 





The measurements of the quantity of air flowing through the col- where h is the pressure drop across tne nozzle measured by the manometer 

d : + : : ` 

lector were performed 5y an accurate technique of measuring the pressure in mm H,O. 3.5 gemenge, bMS = دمه‎ 
drop across a nozzle. The nozzle and the corresponding track shown in Fig. 6 

were fabricated by the firm DEBRO** according to the specifications of 3.4 Radiation Measurenente 








Temperatures are being measured with the help of Pt 100 an 
Pt 500 resistance temperature sensors which have been fixed in a coppe 
with a plug at the other end in a manner shown in Fig. 7a and 7b. As i 
evident a 4 wire technique has been used so that one does not have to 
account for the resistance of the wire leads. The temperature value co 
ponding to the resistance of the sensor is tabulated in the DIN descri 
number 43760 /32/ and follow the following formula: 


The total solar insolation on the collector is measured by an 
Eppley pyranometer (Mode?! PSP) by mounting it at the same inclination as 
the collector. Precautions have been taken to avoid subjecting the instru- 
ment to mechanical shocks or vibrations during the installation. The pyra- 
nometer 1S orientated so that the emerging leads (connector) lie north of 
مدږ مي ېب ټ‎ the receiving surface. The pyranometer is supplied with shields which are 
funnel of tr adjusted parallel to the plane of thermopile and lie just below it. Eppley 


y 
pP ٠ | Pyranometer data are given below. . له‎ YY”. se”) 





Be plexiglass or steel tube P 
/ : | valid in the temperature range from O to 850 9C, where R, is the resis 
/ Instrument characteristics in Mat the temperature ta t is the temperature in 9C and the constant: 





















































Pr 
E 
> 1 = bo are A = 3.90802-1073 deg^!, 8 = 0.580195-1076 deg 2 and Ro = 100%, wh: 
: Sensitivity 10.58 microvolts per watt meter ? is the resistance of the PT 100 sensors at O °C (500 اک‎ for PT 500). 
= Impedance 668 ohms at 23 °C Neglecting the square term in equ. (2) one gets in the 
> linear approximation for the temperature 
Receiver circular | cm ^, coated with Parsons' 
black optical lacquer | l e 
ZER 
O 


Temperature dependence | per cent over ambient temperature range 
-20 to +40 °C (temperature compensation 
of sensitivity can be supplied over other 
ranges at additional charge) 


Using equ. (3) in (2) one gets the expression 


Fig.6: Track specifications for air flow measurements 


Li i m" TD d 
linearity 0.5 per cent from 0 to i400 watts m ` um i "" 8 


+ = 
"e ët 
This rather simple relation is a very good approximation to the exact 1 


E | formula (2), as shown in /33/. The difference evaluated by equ. (4) to 
Deutsche Industrie Norm exact formula equ. (2) is, at t = 65 OC, only At = 1.7۰1073 °C and at 1 











| e i 1 sir heaters st WA-GIE " It is our experience that measuring the air speed with an anemometer 
Fig. 5: A view of the test facility for long plastic solar and then multiplying it by the density of air and cross section of the 
3 outlet to get the rate of the flowing air in kg/s results in overestimates. 
Thus often high efficiencies have been reported in literature /16, 29/. 
DEBRO-MeBtechnik, 4005 Meerbusch !, West Germany _ ! Lx Are | | | | ipe : 
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5 9C): 
ry stan- 
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(4) 


IN- 
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ower 


Multimeter 
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Fig./a: Arrangement for 
temperature sensor (PT 500) 
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— multi-pole plug 
(weather-proof) 


—— mounting screw 


—- metal tube with 4 
connecting wires 


inside 


— PT 500 element with 


protecting wire mesh 


Fig.7b: Photograph of the used PT 500 


assembly 
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temperatures even less. Equ. (4) was used for calculating the temperatures 
of our experiment via the MADAS data processing system (see section 3.8). 


The maximum deviation of the DIN temperature sensors PT 100 
and PT 500 up to 60 °C is + 0.6 °C, at lower temperatures less (At = 
+0.3 + 0.05 [t]). For measuring the temperature in 4-wire technique a con- 
stant current of | mA is sent through the sensor in place (2 wires) and 
the voltage drop across its resistance Rt is measured with the remaining 
two wires giving thus the resistance of Re at temperature t. The change in 
resistance of the PT 100 sensor with the temperature is 0.39 - 0.38.0. /deg 
(approx. 2.5 degrees perf.) in the temperature range 0 °C - 60 OC (approx. 
1.95f1/deg or 0.51 deg/Mfor PT 500). 


| The advantage of using PT 500 instead of the common PT 100 is 
obviously the higher voltage drop (approx. 550 mV) and the higher change in 
voltage due to the temperature change, resulting in a better measuring accu- 
racy. A certain disadvantage may be, that the constant current sources 
(1 mA) must be capable of providing a higher voltage, if a large number of 
PT 500 are used in series. 


- Temperature Measurements in the Collector - 


The inlet and outlet tempratures are measured by fixing the 


sensors at the respective ends of the collector. At the outlet it is expected 


that a temperature stratification occurs. The sensors are therefore so ad- 
Justed as to measure a temperature profile. Five sensors are fixed at the 
outlet while three sensors are fixed at the inlet. An effort was also made 


to mix the air at the outlet by making an arrangement so that the upper layers 


of the air flow to the bottom side and the lower layers flow to the upper 
side. The temperature distribution along the length of the collector was 
also measured by putting five temperature sensors between the inlet and the 
outlet approximately at equal distances from each other. The distribution 
of the temperature was measured above as well as below the absorber sur- 
face (see Fig. 4). 


3.6 Relative Humidity Measurements 





| To get the relative humidity we measured the ambient temperature 
using a dry bulb and also a wet bulb sensor which are then used for the cal- 
culation of the relative humidity from the formula (8). For the ambient 
temperature the dry bulb temperature of the System was taken. The ambient 
temperature sensor is housed in a well ventilated instrumentation shelter 
with its bottom 1.25 m from the terrace and with its door facing north, so 
that the sun's direct beam does not fall upon the sensor when the door 
1s opened. The instrument shelter is painted white from the outside. Precau- 
tion 18 taken not to have any obstructions higher than twice the height of 
the house nearby. 


Relative humidity of the atmospheric air is calculated from the 
measurement of dry and wet bulb temperature using the following relations, 
which are given in /34/. Magnus' formula gives the saturated pressure Po 





* i 
For Tg <0 °C one has to wait with the measurement until the whole wet 


bulb is frozen and equ. (5) has to be replaced by 
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The real vapour pressure is given by Sprung's formula for the ideal psy- 
chrometer as a function of the dry (T4) and wet bulb (T S) temperatures and 


the atmospheric pressure Py 
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where C, is the specific heat for dry air, L, is the heat of evaporation 
for water and P_(T,,) is found by replacing Ta by T,, in equ. (5) (0.622 
is the ratio of the mol-masses of water and air air). 


The constant before the psychrometric temperature difference 
(T, - T,,) in equ. (6) has the value of 0.492 for Pa = 755 Torr and L,, = 
592 cal/g, valid for 10 9C. By convention one takes the constant as 0.5 
without looking at the change of Ly with temperature. 


The relative humidity r as function of the measured dry bulb Tg 
and wet bulb temperature Ty is given as the quotient of equ. (6) and (5). 


P 
z (7) 
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(8) 





Equ. (8) has been used for calculating the instanianeous relative humidity 
of the ambient air blown through the collector. 


The specific heat of the humid air also changes and can be ex- 
pressed according to /35/ as 


Cp (T) - Cpq (T) + W(T)-C,, (T) (9) 


where Cpg is the specific heat for dry air and Con of moisture, which can be 
taken from tables. 


The humidity W in (9) is given by 


(10) 





where r is the relative humidity from equ. (8), P.(T) is the saturated va- 
pour pressure at temperature T and P, is the atmospheric pressure. 


All the quantities in equ. (9) show a small variation with 
temperature in the considered temperature range )10 9C - 65 9C), the effect 
on the performance data of the collector will be discussed in section 6, 


3.7 Wind Speed Measurements 





The wind speed was measured almost at the same altitude as the 
collector test stand with the help of an electronic anemometer obtained from 
the firm Thies /36/. The instrument gives a DC-voltage in the mV- 
range and was calibrated in a wind channel. The experimentally formed rela- 
tion between wind speed v and DC-voltage U is given by the formula /37/ 


v = 0.639 = + 0.04119 وه‎ U (11) 


The smallest voltage difference recorded by the digital voltmeter is 0.1 mV. 


The constant of 0.639 m/s in equ. (11) is due to some turning losses in 
the mechanism of the anemometers. The instrument can be used for measuring 
instantaneous wind speeds in the range of 0.65 - 35 m/sec. 


Data Processing System 





I: every experiment of thermal application of solar energy the 
solar radiation data, temperatures, heat fluxes and other weather data have 
to be measured. These physical quantities are measured by sensors producing 
a voltage or are characterized by a resistance. As the quantities are all 
changing slowly with time (time constants are in the region of a second) one 
can use electronic systems with a low time solution, but having a high 
accuracy converting the signal into digital values. 


The typical voltage values for a solar measuring svstem will be 
U > 100 «V, the time constants of the sensors will be in the range of 
t »100 msec, and the measuring range will be y 1:1000. Then the measurement 
values can be recorded electronically by an analog-digital converter (ADC) 
with a 12 bit solution (1:4096) and a voltage solution of about 100 AV. 
The temperature measurement with a resistance thermometer (PT 100) ‚for in- 
stance, are then performed with an accuracy of abcut 0.2 °K. 


The measurement signals from the different sources are processed 
by the Microprocessor Aided Data Acquisition System (MADAS) designed and 
constructed for the different requirements of meteorological colle ^or and 
system studies /38/. 


The MADAS-System, shown in Fig. 8a + 8b is based on the CAMAC 
standard, which is a convention for a standard communication between differ- 
ent modules and providing a standard power supply to themodules. The main 
component of the system is the MACAMAC (Micro Autonomous Computer Applied 
to Measurement And Controi), an autonomous Crate Controller equipped with à 
8080 INTEL Microprocessor. The controller has an autonomous operation 
system and uses the BASIC language. 


A battery buffered clock is used for controlling the internal 
system functions. Sensors providing analog signals are connected to the 
analog/digital converter (ADC). The ADC has 32 channels, and in build 
scanning system, programmable sensitivity and a constant current source 
for passive sensors. The scanning rate is about 3 sec/32 channels. The 
working of the ADCs is based on the dual slope principal by which noise 
signals are extensively suppressed. Sensors providing pulsed signals (as 
flow meters) can be connected to the counter. Controlled by the clock the 
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Fig. 8a: The MADAS data processing system 


data are shifted to the memory part of the system where they are submitted 
to the combining and averaging procedure before storage in the internal 
memory of the system. For our purpose, the system was used together with a 
HP-Terminal (Hewlett Packard type 2648) where the Basic Programme and the 
data could be displayed. 


In our experiment we used a pyranometer, a wind speed meter 
and various temperature sensors for the dry and wet bulb temperature (rel. 
humidity), for the temperature profile at the inlet and outlet of the collec- 
tor's cross-section, and for the temperature distribution along the length 
of collector. All instruments deliver analogous DC levels which are fed in- 
Lo the ADC. The data recording and storage is controlled by BASIC programme. 
The programme calculates iteratively 10 min mean values and variances of all 
observed quantities which together with time and date are written on a mag- 
netic tape cassette and simultaneously displayed on the terminal. 


The pressure dropAp across the nozzle and the static pressure 
is read from a Betz manometer and fed into the programme via the terminal 
nanually. The programme calculates the air density Ø as a function of the in- 
let air temperature and the static pressure (see equ. (32) and then uses Ap 
ind p for evaluating the air mass flow m (see equ. (1)). From the 1O min 
(or longer) mean values of the inlet and outlet temperature sensors m 
hich is kept constant for longer time intervals and Cp (assumed a constant) 
ind the mean value of the solar radiation the collector efficiency valid 
or the time interval (usually 10 min) is calculated (see equ. (12)). 


| PERFORMANCE EQUATIONS 


Thermai efficiency of a solar collector is usually defined as 





e actual useful energy collected 
1 solar energy incident/intercepted by the collector 
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Fig. 8b: The MADAS-System with plug board and cables from the air collector 
and meteorological instruments. At the bottom is a built in 
Polycomb recorder. 





For each 10 minute segment for which an efficiency value is to be determined, 
the value is calculated using the equation 


ti : 
ise pS, Us Ti) dr 
A. ['r dr 





The quantities m and Cp have been taken out of integration in the numerator 
since they remain essentially constant during one set of measurements. 


Our measurements show (see section 6) that for a constant flow 
rate the variation of average temperature rise, i.e. 


Uu 
1 ١ 
<åT > = €t) Y Í (uch) dt 


with the average incident global radiation, i.e. 


follows approximately a straight line satisfying the equation 


at= 9% / 1% (13) 
nC, 


From equ. (13) it is obvious that the slope of the straight line enables to 
calculate the efficiency. We made a least square fit of the (A'T»vs. <I> 
measurement to get the linear coefficients and hence the efficiency. 


A linear increase of (AT? with <I> for a constant flow-rate 
over a day indicates that the heat loss coefficient from the collectors 
system as a whole remains approximately constant. 


Hottel-Whillier-Bliss equation /39, 40/ states the instantaneous 
performance for an element of solar collector. The equation states that 
the useful energy collected per unit area of the collector is the difference 
between the amount of solar energy absorbed and the heat loss to the sur- 
roundings. Mathematically 


Å .p(eeL- WITT) 


Using the definition of the efficiency one can write 


Up 
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For obtaining the detailed information about the performance of 
collectors and to prevent the necessity of determining some average plate 
temperature, it is convenient to introduce a parameter F' as 


actual useful energy collected 
useful energy collected if the entire surface collector 
were at the average fluid temperature 





F' = 


Mathematically 


du 


T ATA) I - U, (T,-T )] 





Introducing F' in equ. (15) one gets 


U 
pv K : D (Te - rj] (17) 


For the mean fluid temperature it is customary to write (T,+tTi)/2; if one 
plots d against 


1 3o *i 


TSG 


a straight line should result. The slope of the straight line is some function 


of Ur, and the intercept on the efficiency axis is some other function of 
(t44)e, one can therefore determine the collectors characteristic parameters, 
i.e. the heat loss coefficient UL and the collector efficiency factor PI, 

It should, however, be noted that an approximation of putting 
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is not essentially a good one, especially for a porous absorber air heater. 
This is verified by our measurements on the three collectors (see section 6). 


Another form of collector performance equation is generally used 
/6,8,9/. This form is derived from the definition of the heat removal factor 
/4/ viz 


F = useful energy collected 
R useful energy collected if the absorber plate were at the 
inlet fluid temperature 





Mathematically 


du 


d de A[(taa) I - D (T;-T D 








20 


Substituting in equ. (14) and rearranging | 
=F, (Ica), - U, (T; - T0] 


For solar air heaters with mat 
form is definitely valid, but does not lea 
plots % against (Tj-T,)/I. This is due to 
for such porous plate absorbers.” 


5 VARIATION OF AIR DENSITY WITH | 


The density of air changes due 
Following Hammecke et al. /42/ we present | 
accurate formula for determining the air d 
temperature. The compressibility Z (p,T) o 


by 
e El 
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Comparing the density P(p,T) with the dens 
and pressure one gets 
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The form of Hottel Bliss equation for ma 
discussed in a later communication. 


In this section the temperatures are def 
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(24), therefore, becomes 


* 
puta (27) 


From Reference /43/ one can take Z as a function of pressure and temperature 
and calculate a from the equation (26). Figures 9a and 9b show (-a) as a 
function of pressure and temperature. The figures show that a maximum value 
of a is obtained at a temperature of 650 °K. In Fig. 10 we show F(p,T) = 
a(p,6509K) -a(p,T) as a function of (650-T). The dependence comes out to be 
linear. One can write, therefore, for a 


a= ESP [1- p ER 


From Fig. 9a one sees that pressure dependence can be represented by B 
Factor A, can be determined by putting T = 650, so that equ. (28) gives 


Ao = a/ (p/p) 


For the pressure of p - 40 atm Ay = — 3.88 x 1074 


The temperature exponent C is the slope d 1n F(p,T)/d In ]650-T]. One gets 
C = 3.4. The factor D comes out to be 0.575. Hence 


e e و‎ 7.4 
a = -3.88 x 10 4 (E. |] = 0,575 pa | ] 
Po 0 


Substitution of equ. (24) into equ. (26) one gets 
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For T = 273 °K and د‎ = 1 - , 

and p. atm, one can take Po 1.29304 kg/m” and Zo 
0.99921 from the tables in /43/. Hence from equ. (30) and (31) one finally 
gets the following equation for the air density 


سان ا M‏ 


P = 1.2923 P p! - 3.88 x 10. 


— Y d 
Po T 


E> ^ a n ef 


3.4 


a - 0.575 [PRT ار‎ 


0 





A 
= 














then multiplying it by the density of air and cross section of the 
outlet to get the rate of the flowing air in kg/s results in overestimates. 
Thus often high efficiencies have been reported in literature /16, 29/. 
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6 RESULTS AND DISCUSSION 


6.1 10 m Collector a 








We first put the 10 m collector a for the test. From 28th Ju 
1981 to 24th August 1981 the measurements of the solar intensity were pe 
formed on a strip chart recorder and the temperatures were recorded manu 
with the help of a resistance measuring meter. From 25th August 1981 the 
data processing system was available for the measurements. The day long ı 
surements of the relevant parameters were performed for various mass flo 
rates. In Fig. 11 we show one such measurement of various temperatures fi 
a mass flow rate of 371.8 m3/h. Fig. 12 shows the corresponding variatioı 
of the instantaneous efficiency (calculated using equ. (12) over 10 minu 
interval). The fluctuations in the outlet temperature, temperature diffe: 
and the efficiency correspond to the variations of the solar intensity. | 
the peaks in the outlet temperature exactly correspond to the peaks in 
solar intensity curve, the peaks in efficiency are slightly shifted. The 
reason for this being a certain time constant for the collector. The dec: 
in temperature of the absorber plate (and hence the outlet temperature) 
takes place at a slower rate than the decrease in the insolation. A plot 
the temperature rise as a function of solar radiation for various air fk 
rates show that the measurements are scattered around a straight line as 


Equ. (32) was used to calculate the variation in air density in the formula 
for the output power and efficiency equ. (12) of the solar collectors. The 
working range is approximately from p = 0.9 atm to 1.11 atm for the static 
pressure and T = 280 °K to 340 °K for the air temperatures. Some typical 
values of air density at various pressures and temperatures are given in 
table 2. 
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Fig.10: F(p,T) as a function of (650-T) and p as parameter 


Table 2: Variation of air density with temperature and pressure (equ. (32)) in shown in Fig. 13. We made a least square fit of the measurements corres- 
f | the working range 
3 —_ ___-—-—— -= 
2 | Density e (p,T ) in kg/m „| 18 m COLLECTOR a 1 
| i e | 1 . ul | " | dote 200 1981 i - 
R 2 : t bat 35 $ flowrate 371,8 ۸ 3 
2 10 2 » 10 Ze | i 1.2832 | 1 i i 1 e d «uen o 65 Lei n 
| 1. 1.2603] 1. 1. i 1. pra oe eg 
| Ba. | eg 
Fig. 9a: Variation of function a with pressure and temperature as parameter 1. 1 2384 | i a 1 K - ost E 1 | ^ 
a. 1.2168 E. i. i. 1. zZ outlet temp. ` T int 4 450 pa 
1 1.1961 | 4. 1 i. : o 2| ee 2 
| i. 1.1761 | 1. i. i. i. ee Ea و‎ > 
40 atm X 1.1568 | 1. i. B. i. = —— 1 

d 1.1382 i. 1. i. ud aa up يی کن‎ n 
: 1.1201 | 1. 1 1. | Hw j loo. 
U í 1 ; 1025 1 : Í 1 | temp. difference ml! ` - d 
0. 1.0855 | 4. i. i. ا‎ ' 

0 1 B. » 1 

Û. 1 E 1 i 

0 1 1. 1 P 

0. i i. i. 1 

0. 1 P 1 B 


خم خد مو حا pa. de‏ حم Si e‏ 


. 4079 








Fig.11: Time variations of various temperatures and solar intentsiy 
on a typical day 
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Fig. 9b: Variation of function a with temperature and pressure as parameter 
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Fig.12: Time variation of the temperature rise, efficiency and solar 
intensity on typical day. 
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Fig.13: Variation of temperature rise with solar intensity for various air 
flow rates 
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ponding to the appropriate line and back calculated the efficiency from the 
slope of the appropriate curve (eqn.13). The calculated efficiencies for various 
flow rates are also tabulated in the same figure 13. Such a curve, as one can 
see, carries maximum information as it gives for a particular rate of air 
flow and for a particular intensity, the temperature difference and the 
efficiency of the collector. The efficiencies may, however, differ from the 
instantaneous efficiency. In Fig. 14 we have plotted the measurements, 

over 10 minutes interval, of efficiency and temperature rise performed in 
almost constant radiation conditions. A dip in the efficiency around 

flow rates of 660 m^/h and 700 m3/h is due to change in the meterological 
conditions and the related thermal losses. The efficiency and pe temperature 
rise curve show a saturating trend beyond an air flow of 600 m-/h. A maximum 
efficiency of around 50 % shown in Fig. 1h is five percent smaller than 
obtained from the least square fit of Fig. 13. The later, however, contains 
errors of averaging etc., however the straight line graphs are found to be 
good from the point of view of comparing the performance of collectors. 
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Fig.14: Variation of efficiency and temperature rise with air flow rate 


6.2 10 m Collector b 





This collector was tested for a period of one week starting from 
2nd September 1981. Each day the flow rate was changed and the measurements 
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Fig.15: Time variation of solar intensity and various temperatures 
on a clear day 


were recorded. A typical curve showing the ten minute variation of outlet, 
inlet and ambient temperatures is given in Fig. 15. In Fig. 16 the corres- 
ponding variation of the temperature rise and the efficiency are shown. On 
this day the ambient conditions were quite good, and the solar intenisty 
was also constant over long periods. The flow rate of 757 m3/h was one of 
the highest one could achieve from the present blower and hence an effi- 
ciency of about 32 % obtained for solar radiation of about 800 w/m? is 
almost the maximum one can obtain from this collector. This efficiency is 
markedly lower than that obtained from the 10 m collector. The reasons for 
getting this low efficiency are discussed later in section 6.4. 


Measurements of AT, the temperature increase for various air 
flow rates and for various magnitudes of solar intensity are shown in 
Fig. 17. The measurements were fitted to a straight line by the least square 
method and from the corresponding slope of the straight line the efficien- 
cies were calculated which are given in Fig. 17. The efficiencies obtained 
by this method are found to differ with the values of the instantaneous 
efficiencies measured over ten minutes interval for approximately constant 
solar radiation. The instantaneous efficiencies are plotted in Fig. 18. 
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Fig.16: Time variation of solar intensity S, temperature rise AT and 
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Fig.17: Variation of temperature rise with solar intensity for various 
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ig.18: Variation of efficiency and temperature rise with air flow rate 


3 20 m Collector 





This collector was also tested for a period of one week start- 
ng from 9th September 1981. The measurements corresponding to a flow rate 
f 697 m3/h and for weather conditions of 10th September 1981 are shown in 
igs. 19 and 20 for the temperature and efficiency respectively. Though a 
emperature rise of about 17 OC is reached in this case, corresponding to 

solar intensity of about 850 w/m?, an efficiency of only about twentynine 
ercent is obtained for this collector. 


'asurements of AT as a function of I over the entire days for various flow 


ites are plotted in Fig. 21 as a function of the solar intensity. The best 
(t lines obtained through a least square fit of the data, corresponding 


the appropriate flow rates are also drawn in the figure and the effi- 
encies for varıous flow rates, as given by the slope of the appropriate 
raight line are also tabulated. Fig. 22 shows the instantaneous efficien- 


measured over a small duration of about one hour for nearly constant 
tensity of ~ 850 w/m2. 
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Fig.21: Variation of temperature rise with solar intensity 


Fig.19: Time variation cf various temperatures and solar intensity 
on a typical day 
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Fig.22: Variation of efficiency and temperature rise with air flow rate 
on a typical day 





inpet: MEA TEID CED NE lap Cat I oa arn A MIS APP Urt PAP Seh SOPE ATOR, د‎ are ipae a tnn LEUTEN N e .۰ه‎ ON DEE ESIRAN ب“ جم > رس‎ T NUS 


Temperature Difference 


6.4 Comparison of the Collectors 





A comparison of the temperatures of the flowing air at the out- 
let, temperature rise obtained in the collectors for various flow rates and 
the instantaneous (equ. (12)) efficiency of the three collectors for various 
air flow rates is tabulated in table 3. It is easily seen from the table that 
the performance of 10 m collector a is much better than the performances of 
either 10 m collector b or the 20 m collector. The efficiency of 10 m collec- 
tor a is almost double than that of 10 m collector b. The efficiency of 20 m 
collector is still low. Even the temperature rise obtained from the IO m 
collector a is more than that obtained by 20 m collector up to an air flow 
rate of 410 m3/h after which it is almost equal to the temperature rise ob- 
tained in the 20 m collector. The efficiency of 20 m collector is seen to be 
less than the efficiency, of 10 m collector b having the same textile. This 
18 because the temperature rise is not in proportion to the increase in area. 
Due to increase in the area of the collector, thermal losses also become 
more (radiative losses not being linear). It is also noted that the effi- 
ciencies calculated by least square fitting the AT versus I (Figs. 13, 17, 21) 
also differ in the same proportion for these collectors. Increasing flow rates 
results in increasing fan power as evident from the last but one column 
of table 3. This consumed electrical power should be subtracted from the 
thermal power q, = n-I-A to get the useful power. This is also given in 
table 3. It should be noted that though qu goes on increasing with the flow 
rate in an approximately exponential manner, the useful power Q is maximum 
at a particular flow rate. e 


The reason for the difference in the performances of various 
collectors is essentially due to two reasons: 


(1) Due to the difference in the absorptance of the black ab- 
sorber tissue, 


(2) Due to the different density of pores in the two absorbers. 


Table 3: Comparison of the thermal performance of the three collectors 
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130 |20.3120.1/20.2136.6/29.9/37.0/16.3| 9.8/16.8 (48. 3/26.5| 22. 7 13966 |2383 14086] 363 |3603|2020|3723 
20.4 21.7 [19.9137.3/31.8|37.4|16.9 |10. 117.5 |48.0|26.2]22.6 |3942 |2356 |4082| 326 |3616 12030 |3756 
20.8 22.6 |18.7|38.2|33.1/36.9|17.4 |10.518.2 |46.625.7|22.2 |3826 12310 |4001| 299 [3527 |2011 |3702 
21.2 |21.6 |19.5/39.4/32.6|38.2]18.2|11.0]18.7 |45.1/24.8|21.1|3698|2237 13798) 245 |3453| 1992/3553 
21.8 121.7 |[18.7|41.233.3]38.8 |19.4 11.620. 1 40.9 |22.3| 19.3 |3557 |2009 |3481| 220 [3137 | 17893261 
21.7 [21.7 |20.0]45.9/33.9]40.6 23.4 |12.2120.6 38.9 18.5] 15.6 13195 |1665 |2813| 164 |3031 | 1501 |2649 
20. 1]21.9]21.9]50.6|39.8]48.0|30.5 |17.9]26.133.4|17.8] 13.0 |2742 |1602 |2345| 112 1263011490 12233 
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The absorptance of the textile 
lector a is fifteen percent more than the 
(Fig. 3b) used in the collector 10 m collec 
crease in the absorptance is, however, not 
difference in the efficiencies of the colle 
Fig. 23 we have plotted the temperature prc 
of the collectors, above as well as below t 
black textile, the temperature is seen to b 
length in the 10 collector b and the 20 mc 
a the temperature of the air above the abso 
(except at very small rates of air flow). T 
collector a from the top are therefore much 
tors. Moreover, Fig. 3b shows that this tex 
whose surface photograph is given in Fig. 3 
in 20 m collector all the air is not able t 
bvilds a pressure at the outlet end of the 
out from a study of the temperature variati 
collector) below the textile absorber in bo 
20 m collector. The temperature is seen to 
till the outlet end (for 10 m collector b) 
to increased air flow from the top of the a 
20 m collector the maximum of the temperatu 
a distance of 15 m from the inlet and then 
the temperature variation (along length) sh 
the reason being the nonuniform porousity o 
from Fig. 3b. A better absorption of solar 
of the textile absorber of 10 m collector a 
good thermal performance. 
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absorber (Fig. 3a) used in 10 m col- | ; " > 
absorptance of the textile absorber | | 
tor b or 20 m collector. This in- 7 

able to explain rather big 

ctors with the two textiles. In 6.5 
files measured along the length 

he textile absorbers. Above the 


Effect of Back Insulation 


The insulation at the back of the collector improves the effi- 


e increasing steadily along the ciency of the collector especially that of 10 m collector a. This effect é 
ollector. In the 10 m collector was studied in two steps. First we put a simple aluminium foil above the X 
rber is nearly equal to ambient 3 cm thick polyethylene insulation (which was at the back of the collector | | : p 
he thermal losses in the 10 m in all the previously described measurements) and measured the average | m rg | | 4 | | Kë 
less than the other two collec- | 
tile is denser than the textile E : 
a. In 10 m collector b and also i -— 
o go through the textile and Table 4: Performance of 10 m collector a with and without aluminium foil 
collector. This is also born over the back insulation for incident solar intensity of 700 Wn’. 
on (along the length of the | E 
th the 10 m collector and the | 3 T 
be increasing with the length | . 
and suddenly dropping down due | | . er E He 
bsorber to the outlet. In the i | p s E i 
re is seen to be reached around : p | | y 
decreasing. In the 10 m collector a, * 
ows a fluctuating behaviour; | "ug | '. | 
the porous absorber as evident We 1 no -. 
radiation and better porosity | á 
are thus responsible for its | ° . : 
20m COLLECTOR 
eMe Mull i e 
| d 
I Polyethylene as back insulation | | eee ti 
II Polyethylene with an aluminium foil as back insulation | AR | . | SE > | y 
هه‎ ` | 
= 750,3 ; 1=676,8 temperature rise of the air flowing at different rates over ten minutes | | د‎ | ei 
States مس‎ | interval. Typical measurements along with the calculated efficiencies are P | | a | 4 
Summarised in table 4 along with the corresponding measurements performed pony x ep . E | E 


without the aluminium foil. For all the flow rates the efficiencies recor- ! 

ded with an aluminium foil show a remarkable improvement. The unexpected 1 ( + 

result of relatively less increase in the efficiency for the lowest flow- 

rate, i.e. 372 m/h, in comparison to higher flowrates is due to the in- |! 1008 e 
creasing thermal losses through the cover for smaller flow rates. In fact HE des Eo 
for still smaller flow rate (~ 270 m3/h) the efficiency observed with alu- ٢ [S | 

minium foil was even lower than the efficiency without aluminium foil due SE | da X 
to a much colder day in the former case. For very small flow rates the air Ze | rre EE 
keeps mainly on the upperside of the absorber and is not able to flow | i 8 ; * | 
through the pores effectively. We notice not much improvement in the effi- 
ciency of 10 m collector b under similar conditions due to the same reason 
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i Though Fig. 2 shows a rather low transmittance of the PVC foil 
in the temperature of interest, yet to make sure that the zuprovement in 

the thermal performance of the 10 m collector a is not due to the reflec- 
tions of IR radiation, which might come through the back PVC foil, we placed 
a 6 cm thick polyethylene insulation below the collector and carried out 

the measurements which are given in table 5. A remarkable result of obtain- 
ing as high an efficiency as 70.6 percent is achieved for an air flowrate 

of about 800 m3/h. Considering the simplicity and low cost of the collector 
the result is really very encouraging. From table 5 also one notices a 
comparatively lower efficiency for the lowest flow rate. The reason for 

this again being the same as discussed earlier, i.e. sticking of air over 
the upper surface of the absorber for low flowrate and hence increased thermal 


losses from the top. It is, therefore, ovious that for such type of collector 
a threshold or minimum airflow rate 3 500 m3/h is desirable. 


Table 5: Performance of 10 m collector a with different back insulations 


(figures in brackets indicate the incident solar radiation in W/m 
on the collector surface) 





ó Temperature Rise °C Efficiency 2 





(m/h) | I o |I 1 II 
420 23.4 | 23.9 | 22.8 0.4905 | 0.4781 
(900) | (748) | (732) 
19.4 | 18.0 | 20.4 0.5257 | 0.5447 
(836) | (682) | (746) 
18.2 | 18.4 | 19.2 | 0.4898 | 0.5421 | 0.6090 
(876) | (800) | (743) 
17.4 | 13.9 | 18.8 | 0.5029 | 0.5560 | 0.6561 
(908) | (656) | (752) 
15.9 | 13.9 | 18.3 | 0.5233 | 0.5588 | 0.6785 
(855) | (700) | (759) 



































3 cm thick polyethylene 
3 cm thick polyethylene with aluminium foil above it 
6 cm thick polyethylene 


6.6 Pressure Drop in the Collector 





A good porous construction of the absorber plate also results 
in a smaller pressure drop across the collector. The pressure drop, measured 
by the static pressure difference across the inlet and the outlet, in the 
10 m collector a and 10 m collector b as a function of various air flow rates 
is shown in Fig. 24. As expected the pressure drop across the 10 m collector b, 


which has the denser absorber textile of Fig. 3b, is much higher 


(approximately a factor of 2) than the pressure drop across the 10m collector a. 
In collector a the air is steadily passing through the absorber along the 

whole length, while in the collector b more air flows along the top up to 

the end and then passes through the textile to the outlet. Since the air 

in textile of Fig. 3a passes chrough it, the heat transfer is from the en- 

tire volume resulting in high average AT and hence the higher efficiency in 


comparison to the corresponding performance of the textile of Fig. 3b. 
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Fig.24: Pressure loss as a function of air flow rate in the two 
collectors of the same length 
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6.7 Effect of Relative Humidity on the Thermal Efficiency 





Relative humidity of the inlet air affects the thermal efficien- 
cy of the solar air heaters through its specific heat, which varies with the 
amount of water vapour present in the air according to the formula given in 
equ. (9). Calculated values of C at various temperatures and relative hu- 
midity are given in table 6; the saturated vapour pressure at various tempera- 
tures has been taken from the standard tables /45/, The specific heats Cpg and CPm 


6: Variation in the specific heat of air with relative humidity r 
and temperature T 
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of dry air and moisture respectively also change with temperature, bui 
variation being slow; viz 


1.006 <c, (T) ک‎ 1011 kJ/kg °K 
d 
in the temperature range 0 T « T « 80 "th 


1860 < c, (T) € 1.884 kJ/kg °K 
m 


in the temperature range 7 a < T < 87 "e 


we have used for the specific heat calculations, the average values 


E 
C 
Pm 


In a small temperature range, P CT) can be expressed by tt 
following linear relationship: 


1.008 kJ/kg °K 


1.871 kJ/kg °K 


p, (T) = Ry + R,T 


Consequently the humidity ratio defined by equ. (10) can be written a: 


R, + R, T 


W = 0.622 r — ag 





P] 


In the temperature range 10 ?C to 70 "b. the values of R, and R5, obt: 
through a least square fit, are given in table 7 along with the corre] 
with the exact values. It is obvious that the maximum change in the v: 
of c, in the temperature and relative humidity range of interest coulc 
about seven percent maximur. The corresponding change in the thermal € 
ciency is also then about the same, i.e. a maximum of seven percent. ] 
also to be noted that these changes are the maximum possible one. Acti 
changes in the efficiency values are still of smaller magnitude. It is 
therefore not essential to incorporate the change of cp with relative 
dity for evaluating the efficiency of the collector, in particular whe 


aim is only to compare the performance of one collector from that of t 
other. 


Table 7: Constants Rı and R2 obtained through a linear lea 
fit of the saturated pressure of water in a certain 
temperature range 


E "Y R, atm R, atm/deg Correlation 
100-300 -4.788-10 ? 1.487.10 ? 0.9888 
309-50? ° -8.316.10 ^ ۹57 0.9920 
509-709 0 -3.616-10 ! 9.484-10 2 0.9930 

















CONCLUSIONS 
We, however, observe that for 10 m collector a, 















































T af, | | 
dens 1 vs )ال‎ > E14 J and vs 1 (a, - T )] plots cannot be approxi- | Flexible mattress type of plastic solar air heater having a 
6.8 Validity of Hottel-Bliss Equation i di سا‎ i Tem | black porous polyester textile yield high thermal efficiency with reasonable 
Zoch ٩ | | . y : : i ae by a 1177 line. OH ` e freie? purposes we give a typical plot in Fig. 25 high temperature increase. The important thing for such air heaters is an op- 
Equs. (17) and (19) are usually employed for characterising the 9/ ۱ | 10m COLLECTOR b Pe e P ue SE 323: m^/h through the collector. For higher flow rates of air timum pore density allowing a good flow of air and a good absorptivity of the 
performance of solar water /35, 44/ and solar air heating collectors /3,6,8, . | also, similar scattering of points is obtained. The reason for this wide scatter of solar radiation. The collector should be well insulated at the back. 
For flat plate solar water heaters and for nonporous absorber type of air Y flow rate = 757 m3/h points is the complicated nature of volumetric heat transfer in the matrix air heaters 
heaters, the factor PI, F, and heat loss coefficient Uj are nearly constant ' FIUL = C,0168 | سر دا‎ gs go: of F' and ۴ on the extinction coefficient of matrix, its thickness, Since the efficiency as high as about seventy percent is ob- 
and a linear representationof Fp’ f = 0,915 a vhs O the absorber material etc. The mean fluid temperature also cannot | tained, collecting nearly .65 kW/m2 for a solar radiation of 900 W/m?, 
T_+T; -] > e taken to be simply equal to with a temperature rise of nearly 20 °C, the collector can be put to a 
^ vs D, WS > r.l I c 8: variety of applications. Apart from crop drying which is a natural appli- 
T llect au 2 cation of any solar air heater, this collector can be used for domestic 
and vs (T,-T,) plots provide a very useful method for nn en A = Even for nonporous absorbers, if they are long solar collectors, the approxi- | hot water heating, space heating and air conditioning. For the later appli- 
parameters like heat loss coefficient, BESET STS np »ايم‎ meet m ها‎ mation Tf cations it is desirable to have a closed loop cycle for the flowing air. 
product, heat flow factor etc. For matrix air heaters also various authors uJ T +r. 
/6,9/ have assumed the validity of eqs. (17) and (19) for characterizing 0,1 T; = 0 - 1 
their collectors. 
à is not strictly valid. The modified form of F', Fg and Tf for porous and 
T i » nonporous solar air heaters shall be described elsewhere. Acknowledgements 
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